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1. Introduction 

Like in the well-studied Drell-Yan process, the partonic cross sections qq,qg,gg — > HH' +X of 
heavy-particle pair production contain higher-order terms 

[a s \n 2 P]", p 2 = \- z = \-M 2 HH ,/s (1.1) 

in their perturbative expansion, which should be summed to all orders, if it can be argued that the 
hadronic cross section is dominated numerically by these threshold logarithms. For two coloured 
particles in the final state with some fixed invariant mass Mhh' complications arise relative to the 
Drell-Yan process due to colour-exchange and kinematics-dependent anomalous dimensions sim- 
ilar to di-jet production [jl]]. For the total partonic cross section the only parametrically enhanced 
logarithms arise from the true production threshold, M HH < = Mh + M H i and the kinematical de- 
pendence disappears. On the other hand, the particles are non-relativistic in this region and even 
more strongly enhanced terms (a s /P) n appear due to the Coulomb force. Although resummation 
for total partonic cross sections has been performed in the past [Q], the issue of factorization of 
soft and Coulomb gluons, and their simultaneous resummation has not been addressed with rigour 
until recently. In this proceedings article we discuss the factorization formula applicable to this 
situation, the diagonal colour basis for the leading soft function relevant to the heavy-particle pair 
production process, and the two-loop anomalous dimension for soft radiation. For details we refer 
to [|3|]. See the talk by P. Falgari [Q] for a discussion of resummation of the squark anti-squark 
production cross section at the LHC. 

To define the LL, NLL, etc. approximations of the resummed cross section in the presence of 
Coulomb effects, we note that near threshold the usual expansion, where a v lnj8 in the exponent 



of (L2) below counts as order one, is combined with an expansion in j3, such that a s / j8 also counts 
as order one. This leads to a parametric representation of the expansion of the cross section in the 
form 



Oppi 



&{ ) L (it) exp lnj3go(«,ln j 8)+g 1 (a $ ln j 8) + a,.g 2 (a,ln j 8) + .. 

(LL) (NLL) (NNLL) 

x { 1 (LL,NLL); a s ,f5 (NNLL); aj, a s f5 ,f5 2 (NNNLL); . . . } , (1.2) 

which reproduces the standard structure ^ away from threshold for k = and no expansion in 
j3. This implies that at 0((X 2 ) relative to the Born cross section the terms a 2 x { 1//3 ;ln 21 j8;/3 x 
In 4,3 j8} are NNLL. Note that a 2 In 2,1 j8 terms may arise from the product of a Coulomb-enhanced 
one-loop correction a s /j8 and a /3-suppressed soft emission a $ j3ln 2,1 j8. This highlights the subtle 
point that one must consider approximations one order beyond the standard eikonal approximation 
to capture all NNLL terms. 

2. Factorization of soft and Coulomb gluons 

The factorization of soft and Coulomb gluons is a non-trivial issue, since the non-relativistic energy 
of the heavy particles is of the same order as the soft gluon momentum. Hence, Coulomb exchanges 
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are not part of the hard process and take place (diagrammatically) "in between" the soft gluon 
emissions. In [J3j] it is shown that both effects are simultaneously resummed by means of the 
formula 

d(P,^) = ZL H ii'(M^) fdco^JKE-^W^ico^), (2.1) 

a i,i> J R a L 

which contains a multiplicative short-distance coefficient H~, in each colour (in higher orders also 
spin) configuration, and a convolution of soft functions W a ,' a with Coulomb functions 7| . The 
factorization of soft gluons from collinear fields follows from a field redefinition with a light-like 
Wilson line in soft-collinear effective theory (SCET) in the usual way ||5|]. To prove the decoupling 
from non-relativistic fields we redefine \\f a (x) 



Sv (x°) a b Yb 1 ' ( x ) w i tn a time-like Wilson line 
dt vA c s (x + vt)T^ c 



J® < 



si R) (x) = Pexp 



-Igs 



This has the effect of turning D° into d° in the leading-order PNRQCD Lagrangian 



(2.2) 



PNRQCD 



\ 2niH 2 1 \ 2m/// 2 



d 3 r 




¥' 



y/^T^'Vl (0). 



(2.3) 



The key observation is that 5 V drops out from the Coulomb interaction expressed in terms of the 
new fields, since in any representation R and since 5 a d in the adjoint 

representation is real and independent of r. For the latter point it is important that the soft gluon field 
in depends only on x° [Q], while the Coulomb interaction is non-local but instantaneous. This 
proves decoupling of soft gluon and Coulomb resummation, since soft gluons disappear from the 
leading-order Lagrangians for the other fields. They do not disappear from higher-order terms in the 
SCET and PNRQCD Lagrangians, but these sub-leading interactions can be treated as perturbations 



in j6, resulting in the sum over a in (2. 1 ). 



3. All-order diagonal colour basis for the leading soft function 

The soft function relevant to the leading power in the j8 expansion (a = 1 in (j2~l|)) is given by the 
Fourier transform of W* a (z,n) = ^} c {a}^{2}( z ^) c {b} wnere 

H0|T[S v *^^ (3-D 

and {a} = ai<32«3«4 denotes a colour multi-index for the 2 — > 2 production process. W Ra (z,jx) 
is obtained from (|3Tj) by projecting it with Pfs on an irreducible representation R a in the prod- 
uct R®R' = '£pRp of final state representations, and on the elements cfK of a suitable colour 
basis. The number of basis elements (/, i' = 1 . . . n) is constrained by colour conservation. Fur- 
ther decomposing the product of initial state colour representations into irreducible ones accord- 
ing to r® / = La r a> the number of basis elements equals the number of pairs P; = (r a ,Ra) of 
equivalent representations r a and Rp. For example, in case ofa8(£>8— ^8(£>8 process, Pj € 
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{(1,1), (8 5 ,8 S ), (8 A ,8 5 ), (8 A ,8 A ), (8 5 ,8 A ), (10,10), (10,10), (27,27)}, so the basis is eight- 
dimensional. One then finds [^] that W Ra ((0,n) in ( |2.l| ) is diagonal to all orders in perturbation 
theory in the colour basis 

JO 1 . r r « r*"* mi 

1 * y/dm{r a ) 

constructed from the Clebsch-Gordan coefficients that couple the initial and final colour represen- 
tations to the equivalent pairs Pi. This result follows from colour conservation; the fact that the 
Coulomb interaction is diagonal in the irreducible colour representations; and the ability to com- 
bine the two final-state Wilson lines into a single one using Ca% ia2 S[ R J ibi S[ R J 2bi = SfAc Ra h ^, since 
the Wilson lines for a heavy-particle pair produced directly at threshold carry the same velocity vec- 
tor v. Finally, due to Bose symmetry of the soft function there is no interference of the production 
from a symmetric and antisymmetric colour octet, which excludes the only possible off-diagonal 
terms after applying the arguments above. 



4. Two-loop anomalous dimension for heavy-particle pair production 

Once the leading soft function is diagonal and reduced to the soft function of an effective 2 — > 1 
process with a coloured final-state particle in representation R a , the soft gluon part of the resumma- 
tion can be done in SCET just as for Drell-Yan production [^]. For NNLL resummation, one only 
needs in addition the two-loop anomalous dimension and the one-loop finite term of the heavy- 



particle soft functions. Since ( p.ip is essentially the "square" of the soft function relevant to the 
effective 2 — > 1 amplitude, the two-loop anomalous dimension satisfies Casimir scaling and can be 
extracted from In particular, potential three-particle colour correlations vanish trivially after 
the reduction to the 2 — > 1 process. The result can be converted to the anomalous dimension re- 
quired in the Mellin-space resummation formalism [Q], which requires the one-loop calculation of 
the soft function, and one finds [^] 

,(!)*„ _ _ /460 4n 2 _\ 176 



Hff = -C Ra C A [~ 9 r + 8C3j+ -^CRjpnf, (4.1) 

which has been confirmed independently in [^]. The process-independent ingredients for NNLL 
resummation of threshold logarithms in arbitrary production processes of heavy coloured particles 
at hadron colliders are now in place. The resummation of the squark anti-squark production cross 
section at the LHC at NLL using the formalism outlined above is discussed in [Eh. 



5. Threshold enhancements at NNLO 

When resummation is not required one can use the fixed-order expansion of ( |Q| ) to calculate the 
velocity-enhanced terms at NNLO. However, to obtain all G^ln ' /3 correction terms, one must in- 
clude logarithms from sub-leading heavy-quark potentials and sub-leading soft gluon couplings [BJ. 
These additional terms have been worked out explicitly in [10]. The contributions from sub-leading 
soft gluon couplings could be a source of velocity-enhanced three-particle colour correlations, 
which have been calculated in [ITT]] and found to be non-zero at the amplitude level. The result 
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of [[K]] implies that there are no contributions to the ln/3 terms from such three-particle correla- 
tions in both, the virtual and real corrections to the total cross section. This holds independent of 
particular colour representations. 

Here we provide the velocity-enhanced terms at NNLO for the production of a pair of heavy 
particles with equal mass m in the scattering of massless partons in colour representations r and 
r', respectively, under the only assumption that the Born cross section admits an 5-wave term 
proportional to /3 . The heavy-particle pair is in colour representation R a and a definite spin state 
(singlet or triplet for spin- 1/2 fermions). Denoting by o% the NNLO correction relative to the 
Born cross section, the threshold expansion reads [10] 

n 2D R a f/ ou. , „ J,_2 f 2m P 2 



.(2) _ 4* 4 Dl 



+ ■ 



(C r + Cy 



3/3* ■ p 

- SC Ra - lax - 4Re + 2/3 ln (^) 

+ 128 {C r + C,-') 2 In 4 jS + 64 (C r + C/) { 4 (C r + C, 



21 



+ 2(A, + 4C R Jln 



/2m/3 : 



2) 



2C Ra ln 3 /3 



+ 



16 



-(C r + Q) 2 [72L^ - 288L 8 + 576 - 357T 2 ] + — (C r + Cy) 



18Re[C 



(1) 



+ 18/30 (-Lg + 2) + 36C Ra (-3L 8 + 7) +C A (67 - 3?r 2 ) 
+ 16C ffK (j3 +2C i?a )lln 2 i3 



20mTf 



8Zi-48Zi + 192 



+ |8(C r + C,-0 2 
+ 2(C,- + Cy 

+ 2C*„ \-24Li +U2L; 



35tt 2 \ 70jt 2 

A L 8 - 384 + ^ + 112^3 



16Re [4 1} ] (-L 8 + 2) + /3 ( -8L§ + 32L 8 - 64 + 



117T J 



.'8 /67 

+ ' 3 (y 



40 ' 



it' 



4 Re [4 !) ] - 4 (/3 + 2C* ) (-L 8 + 3) + C A 



-224 + 
4024 



357T 2 
597T 2 



27 



+ — F— + 28^3 + 



AniT f 



296 

-40L 8 + — - if 



98 27T 2 

y + — 



4£ 



+ 167T 2 D Sa C A -2D Sa (l+v : 



spin ; 



lnj8 + ^(i; 



(5.1) 



Here C r , Cy and C Ra denote the quadratic Casimir operators of the colour representations, /3o = 
yC^ — \n\Tf is the one-loop beta-function coefficient, and L 8 = ln(8m/ju). The quantities D Ra , 
a\ = — Y n iTf and v sp i n are connected with the heavy-quark potentials such that v sp i n = and 
—2/3 for a pair of spin- 1/2 fermions in a spin-singlet and spin-triplet state, respectively, and D% a 
refers to the strength of the Coulomb potential in representation R a (Dn a = —Cp for the singlet 
and Dr = — (Cp — C A /2) for the octet representation). The only process-specific input is Re [C%\ 
equal to one half of the the one-loop short-distance coefficient Hx{M,pL) in ( |2.1[ ), when the heavy- 
particle pair is in colour and spin state X. Alternative to a direct computation, it can also be deduced 
from the constant term in the threshold limit of the NLO production cross section in colour 
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and spin channel X by comparing the expansion of 0% ' to the formula 




- 4 (C Ra +4(C r + O)) In ( + \2C Ra + 2Re [C?\ + 



(5.2) 
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